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INTRODUCTION 


Since  this  final  report  has  some  historical  flavor,  the  idea 
of  measuring  the  occupation  of  hydrogen  traps  with  positron 
annihilation  came  during  a  TMS  meeting  on  hydrogen  embrittlement 
and  after  the  Principal  Investigator  had  worked  at  BNL  studying 
vacancies  in  NiO  with  positron  annihilation  -  basically  if  a  known 
trap  for  a  positron  was  occupied  by  hydrogen,  the  positron  would 
no  longer  be  trapped  or  annihilated  in  it  because  of  the 
similarity  of  charge.  Thus  positrons  could  give  valuable 
information  about  traps. 

This  final  report  covers  research  carried  out  first  at 
Northwestern  University  in  the  Materials  Sciences  Department  where 
Yong-Ki  Park  obtained  his  Ph.  D.  He  spent  a  summer  at  Brookhaven 

working  with  Drs.  Kelvin  Lynn  and  Louis  Snead  as  well  as  with 

three  of  their  collaborators.  Prof.  Innes  McKenzie  (Univ.  of 
Guelph),  A.  Vehanen  and  B.  Nielsen  (Ris^  National  Lab,  Denmark). 
Subsequent  positron  research  was  conducted  in  the  Physics 
Department  of  Michigan  Technological  University. 

To  put  things  into  perspective.  Prof.  Doyama  (Tokyo 

University)  and  Cotterill^  stated,  -:^,There  is  not  much  of  an 
opening  of  ions  at  the  core  of  the  dislocations.  Therefore 

positrons  probably  do  not  annihilate  at  the  normal  sites  of  the 
core  of  the  dislocations.-  The  idea  that  there  was  insufficient 
room  for  a  positron  to  be  trapped  at  a  dislocation  was  based  on 
some  approximate  calculations  carried  out  while  they  were 
associated  with  Argonne  National  Lab.,Jn  their  opinion,  positron 
might  be  trapped  at  a  jog  because  "...it  looked  like  a  vacancy..." 
and  vacancy  trapping  was  already  known.  The  latter  in  a-iron  had  a 
lifetime  of  ca.  175  psec.  The  mean  lifetime  which  Yong-Ki  Park 
found  that  summer  was  165  psec  and  in  some  cases,  155  psec.  It  was 
Vehanen’ s  surmise  that  the  value  of  165  came  from  a  mixture  of 
vacancies  and  dislocations  with  the  latter  being  due  to  an 


M.  Doyama  and  R.  M.  J.  Cotterill,  in  Proceedings  of  the  5'"*' 
International  Conference  on  Positron  Annihilation,  Japan 
1979.  eds,  R.  R.  Hasiguti  and  K.  Fujiwara  (Japan  Inst,  of 
Metals,  Sendai)  p.89 
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unspecified  type  of  dislocation. 


Vfhen  Park  returned  from  Brookhaven,  working  with  single 
crystals  supplied  by  Prof.  Meshii  at  Northwestern,  he  found  that 
under  conditions  of  loading  designed  to  produce  screw  dislocations 
in  BCC  iron,  the  lifetime  was  as  low  as  142  psec.  With  further 
careful  work,  he  was  able  to  identify  separate  positron  lifetimes 
of  edge  and  screw  dislocations  in  single  crystals  of  a-iron. 

Prior  to  this  time  and  specifically  at  the  Arlington 

International  Conference  on  Positron  Annihilation  in  1981  thr 

Doyama-Cotterill  opinion  was  still  widely  accepted  and  the 

2 

positron  group  at  Argonne  National  Lab  published  calculations 
dealing  with  the  thermal  de-trapping  of  positrons  from  jogs  in 
FCC  aluminum.  It  was  based  on  the  idea  that  positrons  would  not  be 
localized  but  rapidly  move  through  dislocation  cores  as  in  pipe 
diffusion.  No  experimental  evidence  was  provided. 

In  view  of  the  controversy  and  in  preparation  for 
publication,  the  lifetimes  in  BCC  iron  were  verified  by  Park'*  with 

(i)  50  runs  using  10  different  specimens  gave 
114  +  2  picosecs  for  annihilation  in  the  bulk 

(ii)  40  runs  using  9  different  specimens  gave 

142  +  5  picosecs  for  annihilation  in  screw 
dislocations 

(iii)  65  runs  using  7  different  specimens  gave 
165  +  3  picosecs  for  annihilation  in  edge 
dislocations 

Positron  lifetime  results  were  then  used  to  determine 
dislocation  densities  and  the  results  were  verified  using 
transmission  electron  microscopy  (TEM)  and  etch  pit  techniques  on 
the  same  crystals.  A  few  confirmatory  experiments  were  made  on 


L.  C.  Smedskjaer,  M.  Manninen  and  M.  J.  Fluss,  J.  Phys.  F,  19 
(1980)  2237. 

Yong-Ki  Park,  Thesis,  Northwestern  University,  1986. 
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annealed  and  on  twisted  iron  whiskers  . 

Internal  friction  measurements  (to  show  the  effects  of 
trapped  hydrogen)  were  also  made  on  annealed,  cold-worked,  and 
charged  high-purity  single-crystal  iron  specitens^.  In  fairness, 
a  subsequent  study  by  Hidalgo  et  al.^  using  one  deformed 
polycrystalline  iron  sample  failed  to  separate  several  lifetimes. 
They  strained  aluminum  samples  in  ways  designed  to  produce  one  or 
more  jogs  and  observed  some  differences.  Hence  they  concluded  that 
positrons  were  not  trapped  by  dislocations  but  by  Jogs.  However, 
they  presented  no  further  supporting  evidence.  They  attributed  the 

7 

mobility  to  pipe  diffusion  oif  the  hydrogen  like  species..  Waber 
has  stated  that  various  experts  are  not  aware  of  any  experimental 
evidence  for  pipe  diffusion  of  interstitial  atoms  in  iron  or 
other  BCC  metals.  However,  it  could  be  argued  that  the  present 
results  only  apply  to  BCC  metals  where  pipe  diffusion  was  not 
readily  observed. 

8 

Recently,  Shirai  et  al.  has  measxired  the  lifetimes  of 
positrons  in  vacancies  and  in  various  forms  of  dislocation  in 
carefully  prepared  Al,  Au,  Cu  and  a  Cu  alloy  specimens.  They  find 
distinct  lifetimes  for  each  type  of  defect  and  the  data  fall  on  a 
general  line  when  normalized.  These  authors  take  exception  to  the 
Doyama-Cotterill  position.  They  believe  that  the  positrons  indeed 
are  trapped  in  the  dislocations  and  that  annihilation  occurs  in 
the  dislocation-core  region.  Their  positron  lifetimes  are  roughly 

4 

Niel  Roth  and  James  T.  Waber,  Unpublished  Results, 
Northwestern  University,  1987. 

^  Yong-Ki  Park,  James  T.  Waber,  M.  Meshii.  C.  L.  Snead,  Jr,  and 
C.  G.  Park,  Phys.  Rev.,  34  B  (1986)  823-836. 

^  C.  Hidalgo,  S.  Linderoth,G.  Gonzales-Doncel  and  J.  San  Juan, 
Proc.  8th  Inti.  Conf.  on  Positron  Annihilation,  Delft,  19989. 
(World  Scientific  Publishers,  Singapore,  1990)  pp.  371-373. 

7 

James  T.  Waber,  Private  Communication,  (1988). 

g 

Y.  Shirai,  K.  Matsumoto,  G.  Kawaguchi  and  M  Yamaguchi,  Proc. 
9th  Inti.  Conf,  Positron  Annihilation,  Szombathely  1991.  ed 
Cs.  Szeles,  (Materials  Science  Forum,  Aedermannsdorf , 
Switzerland,  1992)  in  press. 
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proportional  to  the  magnitude  of  the  Burgers’  vector.  So  the 
results  of  Hidalgo  et  al.  are  not  supported. 

While  Hirth  had  estimated  binding  energies  of  hydrogen 

9 

atoms  to  edge  and  to  screw  dislocations  in  his  Campbell  Lecture  , 
the  development  of  a  non-destructive  method  of  determining 
dislocation  concentrations  afforded  a  new  approach  to  identifying 
hydrogen  traps  and  determining  whether  they  were  occupied.  It  also 
became  possible  to  estimate  the  binding  energy  of  a  hydrogen  atom 
in  a  dislocation  from  the  measured  diffusion  coefficient.  Park 
and  Waber^^  showed  that  after  the  specimens  had  been 
electrolytically  charged  with  hydrogen,  the  trapping  of  positrons 
by  edge  dislocations  was  reduced  by  20  percent  and  by  about  40 
percent  in  the  case  of  screw  dislocations.  Trapping  returned 
after  holding  a  specimen  for  several  hours  in  a  vacuum  desiccator. 

The  research  then  developed  along  two  paths:  one,  was  to 
study  more  about  the  behavior  of  dislocations  and  the  other  was  to 
study  the  deposition  and  diffusion  processes  for  hydrogen. 

II.  Further  Positron  Annihilation  Studies 

Waber  and  Jong-Lam  Li,  a  Post-doctoral  Fellow  from  the 
Korean  Advanced  Institute  of  Science  and  Technology,  pursued  this 
non-destructive  method  of  determining  edge  and  screw  dislocation 
densities  at  Michigan  Tech. 

Dislocation  Motion 

One  of  the  important  questions  concerning  dislocations 
is  how  they  are  removed  during  thermal  annealing.  It  is  clear 
that  many  of  them  can  be  annihilated  during  the  recovery  stcige; 
more  complex  processes,  namely  recrystallization  and  grain  growth 
occur  at  later  stages  and  exhibit  further  effects.  The  text  books 
all  illustrate  the  disappearance  of  dislocations  by  the  diffusion 
of  vacancies  to  the  extra  half -plane  on  the  slip  plane  where  the 
former  is  "eaten  away"  by  arriving  vacancies  and  climb  occurs  to 

^  J.  Hirth,  Hetall.  Trans.,  A  11  (1980)  861-869. 

Yong  Ki  Park  with  Prof.  Waber,  See  references  3  and  5. 
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new  slip  planes.  To  "dissolve"  a  half  plane  would  required  vacancy 
concentrations  in  the  vicinity  of  lO"^  molar.  However,  it  was 
observed  that  the  lifetime  of  175  psec  which  is  associated  with 
vacancies  in  a-iron,  was  not  a  significant  component  of  the 
measured  lifetimes  of  cold-worked  samples.  While  vacancies  might 
be  generated  by  dislocation  interactions  during  deformation,  few 
of  them  remained  afterwards.  Thus  it  was  unlikely  that  after  the 
production  of  vacancies  by  the  interacting  dislocations  (the 
former  were  assumed  to  have  diffused  away  to  free  surfaces),  the 
process  could  be  reversed  and  dislocations  could  be  reduced  in 
number.  However,  Park  et  aJ.^^  did  not  detect  a  significant 
number  of  vacancies  either  before  or  after  annealing,  although  the 
method  is  inherently  capable  of  detecting  vacancies  at  this  level. 

The  study  of  this  fundamental  metallurgical  question  was 

undertaken  by  Jong-Lam  Lee  and  Waber  using  these  positron 

annihilation  techniques  and  deformed  single  crystals  to  prove  Li’s 
12 

theory  of  recovery  processes  -  namely  that  a  vacancy  mechanism 
was  not  dominant  in  the  thermal  elimination  of  dislocations  -  but 
that  recovery  occurred  by  the  annihilation  of  dislocation  dipoles 
on  glide  planes.  Two  separate  activation  energies  were  determined 
both  of  which  were  much  smaller  than  that  for  self-diffusion  of 
iron  atoms  -  a  crucial  component  of  the  conventional  vacancy-based 
theory.  Screw  dipoles  had  the  smaller  activation  energy  and  in 
Isochronal  runs,  they  had  virtually  disappeared  (at  lower 
temperatures)  before  significant  reduction  of  edge  dipoles 
occurred. 

13 

Specifically,  Lee,  Waber  and  coworkers  developed  two 
kinetic  expressions  which  relate  the  amount  of  change  during  each 
isochronal  heating  step  to  the  overall  kinetics  and  the  activation 
energies.  To  verify  the  rates,  several  isothermal  runs  were  made 


Park  et  al,  See  References  5  and  3. 

J.  C.  M.  Li,  in  Recrystallization,  Grain  Growth  and  Texture, 
(Amer.  Soc.  Metals,  Metals  Park,  1966)  p.  45ff 

Jong  Lam  Lee,  James  T.  Waber,  Yong-Ki  Park  and  Jeff  M.  Th. 
de  Hosson,  Mater.  Sci.  Eng.,  81  (1986)  379-390 
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with  the  same  crystals  and  rates  were  in  good  agreement  with 
the  isochronal  data.  These  isothermal  runs  were  important  in 
establishing  that  the  process  was  governed  by  a  second  order 
kinetic  equation  i.e.,  depended  on  the  product  of  p* p~  and  not  on 
one  of  them  alone,  as  would  be  the  case  with  the  vacancy  diffusion 
mechanism. 


Two  major  conclusions  of  this  study  were:  (1)  that  at 
low  strains,  the  dislocations  exist  primarily  as  dipoles  and  glide 
on  common  planes  with  low  activation  energy;  the  screw  components 
are  more  mobile  and  have  the  lower  activation  energy,  auid  (2)  the 
excess  density  which  is  equal  to  Ip*-  p~i  is  a  small  fraction 
of  the  total  density  (p*  +  p  )  or  p  As  predicted  by  Li,  the 

O 

concentration  of  screw  dislocations  must  become  small  before  the 
edge  components  decrease  significantly. 

14 

After  reading  a  paper  by  Marcinkowski  on  surface  dislocation 
dipoles,  Waber  realized  that  positrons  could  be  used  to  determine 
their  surface  concentration.  Lee  and  Waber^^  undertook  a  study  of 
the  effect  of  metal lographic  polishing  on  the  distribution  of 
dislocations  in  sub-surface  layers  of  iron  single  crystals  i.e., 
the  depth  profile  of  edge  and  screw  dislocations  produced  by 
mechanical  polishing.  This  work  has  general  Importance  since 
almost  all  of  the  samples  used  in  a  wide  variety  of  experimental 
studies  of  surface  phenomena  are  mechanically  polished. 

It  was  generally  assumed  that  dislocations  would  readily 
diffuse  to  and  escape  at  any  free  surface.  Contrary  to 
expectation,  the  dislocation  densities  were  larger  at  the  specimen 
surface.  With  normal  polishing  the  dislocations  penetrate  into 
the  first  50  microns  of  well  annealed  single  crystals  -  even  using 
0.05  pm  alumina  slurries,  the  depth  is  15-20  microns.  As  shown  in 
Fig.  1,  the  dislocation  density  decreases  logarithmically  with 
depth.  The  thermal  elimination  of  both  edge  and  screw  dislocations 


14  R.  J.  Marcinkowski,  J.  Mater.  Scl.,  19  (1984)  1296  ff 

Jong  Lam  Lee  and  James  T.  Waber,  Metall.  Trans.,  21A  (1990) 
2037-2045 
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Dislocation  Density  (  I/m* 


Fig.  1  Dislocation  Profile  in  Polished 
Single  Crystals  of  BCC  Iron 


It  recently  came  to  the  authors’  attention  that  in  1984 
Kramer,  Feng  and  Arsenault^^  had  determined  : 


(a) -  two  activation  energies  in  connection  with 

disappearance  of  dislocations  in  both  FCC  and 
BCC  metals  and 

(b) -  they  found  that  the  dislocation  density 

decreased  logarithmically  with  depth  from 
the  free  surface  as  their  Fig.  9  shows. 

(c)  They  found  that  the  elimination  of 
dislocations  was  a  two  stage  process.  The 
activation  energy  in  Stage  1  was  14,3 
kJ/mole  and  in  Stage  2,  was  128  kJ/mole.  The 
first  stage  is  rapid.  Stage  2  begins  after 
the  initial  density  p  has  'decreased  about  80 
percent. 

This  was  deduced  using  primarily  Xray  diffraction  line 
broadening  techniques  in  their  study  of  creep.  The  similarity  to 
the  positron  annihilation  study  is  striking.-  the  lower  activation 
energy  associated  with  the  glide  of  screw  components  was  26.6  kJ/ 
mole  and  that  for  edge  was  53.0  kJ/mole.  The  authors’  data  is 
consistent  with  two/third  (or  more)  of  the  initial  total 
dislocation  content  disappearing  before  the  second  stage 
(associated  with  the  edge  dislocations)  becomes  dominant. 

17 

Lee  and  Waber  undertook  a  study  of  the  effect  of 

interstitial  solutes  in  the  single  crystals  on  dislocation 

mobility  using  isochronal  annealing.  Dissolved  hydrogen  increased 

the  mobility  of  the  screw  components  but  delayed  the  movement  of 

the  edge  dislocations  until  higher  temperatures  were  reached.  This 

is  consistent  with  early  work  of  Birnbaum  (Univ.  Illinois).  More 

18 

recently,  he  is  of  the  opinion  that  the  mobility  of  both 

E.  R.  Kramer,  C.  R.  Feng  and  R.  Arsenault,  Hetall.  Trans. 

15A  (1984)  1571-1577.  See  also  Hater.  Sci.  and  Eng.  ,  86 
(1987)  1111-1119. 

Jong-Lam  Lee  and  James  T.  Waber,  Unpublished  Research,  1990 
H.  Birnbaum,  Private  Communication,  (1989' 


16 

17 

18 
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components  Is  comparable.  One  difference  in  procedure,  is  that  his 
observations  were  made  in  the  presence  of  the  electron  beam  of  the 
TEM  microscope,  whereas  in  the  present  work,  the  dissolved 
hydrogen  was  removed  in  vacuo  at  liquid  nitrogen  temperature  after 
each  step  of  isochronal  hydrogen  annealing  but  before  each  set  of 
positron  measurements  was  made. 

Doping  with  amounts  of  nitrogen,  insufficient  to  lead  to 
nitride  precipitation,  reduced  the  glide  of  the  dipoles.  However, 
this  study  was  complicated  by  the  desorption  of  the  nitrogen  from 
the  dislocation  traps  in  the  interesting  temperature  region  and 
is  being  revised.  To  show  the  consistency  of  annealing  behavior, 
one  of  the  expressions  developed  for  isochronal  annealing  is 
plotted  in  Fig.  2  versus  a  modified  time.  Data  for  isothermal  and 
isochronal  are  presented  in  the  left  hand  panel  together  with  the 
behavior  of  the  mechanically  polished  samples.  In  the  right 
panel,  the  same  isothermal  and  isochronal  data  are  plotted  with 
the  data  for  nitrogen  doped  specimens.  These  four  different  types 
of  runs  fall  on  the  same  line. 

Hydrogen  Deposition  and  Diffusion 

In  the  period  of  1970  to  1982,  a  series  of  models  for  the 
diffusion  of  hydrogen  in  iron  or  steel  containing  saturable  and/or 

unsaturable  traps  were  developed.  These  range  from  McNabb  & 

19  20 

Forster  to  lino  Without  going  into  specific  details  of 
individual  papers,  the  over-all  situation  can  be  summarized  as 
follows  : 


(i)  different  assumptions  and  different 
mathematical  techniques  were  used 

(ii)  authors  concentrated  on  simple  parameters 
such  as  time-lag  or  break-through  times 

(iii)  detailed  comparisons  between 

theoretical  expressions  and  experimental 
charging  curve  data  were  scant 


A.  M.  M'^Nabb  and  P.  K. Forster,  Trans.  Met.  Soc.  AIME,  227 
(1963)  618-627. 

M.  lino,  Acta  Met.,  30  (1982)  367-375. 
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CHANGES  IN  DISLOCATION  DENSITY  DISLOCATION  CONTENT  CHANGES  IN  PURE 

FOR  POLISHED  AND  FOR  COLD -ROLLED  SPECIMENS  AND  NITROGEN  DOPED  IRON  SINGLE  CRYSTALS 


Fig.  2  Comparison  of  the  Kinetics  of  Dislocation 


(iv)  little  effort  was  made  to  identify  the  traps 
or  to  quantify  them  -  impure  metal 
and/or  complex  alloy  samples  were  studied 

21 

Researchers  from  New  Zealand  analyzed  a  two-step  process  of 
deposition  at  cathodic  potentials  followed  by  a  rapid  reversal  of 
the  potential  to  "de-plate"  the  adsorbed  hydrogen.  They  considered 
following  model  conditions: 


I  .  Surface  Reaction  Boundary  Conditions  For 

Hydrogen  Deposition 


Reasonable  Models 


(A)  Constant  Surface  Concentration 

Potentiostatic 

(B)  Constant  Surface  Flux 

Galvanostatic 


(C)  Flux  Limited  by  Reactions 
in  Solution 


They  did  some  work  with  deformed  polycrystalline  iron 

22 

specimens  since  they  noted  that  experiments  by  Pumphrey  in  1980 
indicated  that 


(i)  Constant  surface  coverage  of  Hydrogen  on 
steel  was  the  preferable  model  -  but- 

(ii)  Sharp  deviations  from  either  models  (A) 
or  (B)  occurred  after  short  charging  times 

(iii)  Deviations  were  in  the  same  sense  whether 
the  experiments  were  carried  out  in  dilute 
H  SO  or  in  NaOH 

2  4 


Pumphrey  concluded  that  the  effects  were  due  to  the  fact  that 
"...diffusion  of  hydrogen  [atoms!  were  reduced  one  or  two  orders 
of  magnitude  by  dislocations." 

23 

Pound  and  Wright  ,  studied  the  ingress  of  H  atoms  into  plane 


R.  M*^Kibben,  R.  M.  Sharp,  D.  A.  Harrington,  B.  G.  Pound  and 
G.  A.  Wright,  Acta  Met., 35  (1987)  253. 

P.  H.  Pumphrey,  Scripta  Met.,  14  (1980)  695  -791. 

B.  G.  Pound,  G.  A.  Wright  and  R.  M.  Sharp,  ibid  p.  263  ff 
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surfaces  of  zone-refined  iron  by  this  “double  step"  potent iostatic 
nethod.  They  observed  diffusion  control,  rather  than  interface 
control,  and  the  influence  of  eliminating  dislocations  by 

annealing,  as  predicted  by  Pumphrey.  They  found  that  trapping 
rate  constant  k  was  five  times  larger  with  the  cold-rolled 

specimens  than  it  was  with  their  annealed  specimens.  They 
concluded  that  the  traps  were  nearly  irreversible.  They  used  an 
unreasonably  large  trap  diameter  and  found  a  much  lower  number  of 
traps  than  the  present  authors. 

Effect  of  Dislocations  on  the  Diffusion  of 
Hydrogen  in  Iron 

In  fairness,  the  Principal  Investigator  recently 

discovered  that  he  had  failed  to  reference  fully  earlier  work  on 

the  effect  of  hydrogen  in  reducing  positron  annihilation  from 

24 

dislocations  carried  out  by  J.  G.  Byrne  and  his  colleagues. 


Replotting  the  diffusivities  determined  in  Yong-Ki’s 
thesis  for  iron  single  crystals  deformed  in  three  different  ways, 
(Waber  and  colleagues  had  used  col-rolled  single  crystal 
specimens,  single  crystals  specimens  bent  in  a  mandrel 1  at  room 
temperature  and  single  crystal  specimens  pulled  in  tension  at  200° 
K)  -  the  apparent  D^  values  fall  on  a  common  curve,  that  is  as 
shown  in  Fig.  3a.  The  measured  dlffusivity  is  inversely 
proportional  to  the  dislocation  density.  While  that  may  initially 
seem  surprising,  all  the  mathematical  models  predict  that  D 

app 

should  be  related  to  the  density  of  traps  by  a  divisor  similar 
to 


D 

app 


D  /  (  1  +  a  N  ) 

0  T 


If  the  second  term  in  the  parentheses  is  large  compared  to  unity. 


F.  Alex,  T.  D.  Hadnagy,  K.  G.  Lynn  and  J.  G.  Byrne,  "Positron 
annihilation  studies  of  hydrogen  embrittlement"  in  Proc. 
Inti.  Conf.  on  Effect  of  Hydrogen  on  the  Behavior  of 
Materials,  Jackson  Hole  Lodge,  1975,  eds  A.  W.  Thompson  and 
I.  M.  Bernstein.  fAIME,  New  York,  1976)  642-650  In  his 
discussion  of  a  paper  by  R.  Gopahengar,  J.  F.  Wallace  and 
R.  A.  Oriani,  Mater,  Scl.  Eng,,  68  (1984-85)191-196  Byrne 
cites  hydrogen/dislocation  studies  on  nickel  and  copper.  His 
comments  appear  in  Hater.  Scl.  Eng.,  80  (1986)  221. 
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Fiz.  3 


Apparent  H  Diffusivities 
Increased  by  Removal 
of  Surface  Layers 
Damaged  by  Polishing 


the  relation  experimentally  observed  would  be  inverse. 

The  development  of  the  non-destructive  method  of 
determining  dislocation  densities  afforded  the  opportunity  of 
readily  verifying  the  relation  experimentally. 

25 

Lee,  Waber  and  Park  studied  the  rate  of  permeation  in 
mechanically  polished  single  crystals.  It  increased  monotonically 
with  the  removal  of  dislocations  in  the  disturbed  zone. 
Dislocations  in  the  first  50  microns  of  a  specimen  950  microns 
thick  reduced  the  diffusivity  tenfold!  Data  are  presented  in  Fig. 
3b. 

This  is  consistent  with  Pumphrey’ s  conclusion  about  the 
importance  of  dislocations  in  the  surface  layers.  This  is  not  what 
one  would  expect  if  pipe  diffusion  were  dominant.  Note  the 
dislocations  are  traps  not  conduits! 

The  famous  Arrhenius  plot  of  Alefeld  and  Vokl  is  presented 
here  as  Fig.  4.  It  show  that  the  diffusivity  of  hydrogen  in  iron 
as  measured  by  a  number  of  authors  falls  very  close  to  a  single 
line  for  the  data  collected  above  roughly  300°  C.  However,  as  the 
test  temperature  approached  room  temperature,  there  was 
considerable  variability  in  D^.  The  diffusivity  obtained  with  the 
present  mechanically  polished  single  crystal  fell  in  the  middle 
of  the  other  reported  data  for  room  temperature.  When  roughly  the 
first  50  microns  of  the  surface  were  removed  chemically,  the 
resulting  D^  value  increased  by  an  order  of  magnitude  and  fell 
slightly  above  the  Arrhenius  line  labelled  D^.  The  hold-up  of 
hydrogen  isotopes  in  dislocation  traps  is  evidenced  by  a  graph 
in  the  next  section. 


Jong  Lam  Lee,  James  T.  Waber  and  Yong  Ki  Park,  Scripta  Met., 
20  (1989)  823-828. 

J.  G.  Alefeld  and  J.  V5kl,  "  Hydrogen  Diffusion  in  Metals" 
in  Hydrogen  in  Metals,  eds.  A.  S.  Nowick  and  J.  J.  Burton 
(Springer  Verlag,  Berlin,  (1983)  p.  231. 
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Diffusion  eotfficient  0 


toV  1en>pef  aturt 


Fig.  4  Diffusion  coefficient  of  Hydrogen  in  a-Fe  (2). 

■  :  hydrogen  diffusivity  in  iron  before  mechanical  polishing. 
A :  hydrogen  diffusivity  in  iron  after  mechanical  polishing. 


lY.  Rapid  Electrolytic  Charge  &  Discharge 

Mobility  of  Hydrogen  and  Deuterium 

27 

Shi  and  Waber  observed  results  similar  to  those  of 
Pumphrey.-  we  found  that  the  trapping  constant  k  was  increased  by 
the  presence  of  dislocations.  Wright  and  coworkers  used  H^S  as  a 
poison  instead  of  and  both  sets  of  researchers  observed  that 

the  surface  coverage  became  nearly  complete  as  the  overvoltage  was 
sufficiently  increased.  However,  Pound  and  Wright  noted  that  the 
coverage  is  quite  small  in  the  absence  of  a  poison. 

Shi  and  Waber  calculated  the  surface  coverage  and  the  total 
number  of  H  or  D  atoms  left  in  the  solid  after  the  cathodic 
potential  was  changed  to  anodic.  Those  atoms  which  diffused  out 
were  oxidized  to  hydrogen  ions.  Then  the  number  of  hydrogens 
remaining  was  related  to  the  surface  density  of  dislocations. 
Half  of  the  residual  atoms  were  assumed  were  to  be  trapped  in  the 
dislocations  and  the  remainder  were  assumed  to  be  in  the  bulk  in 
interstitial  positions.  -  this  is  a  conservative  estimate  in  view 
of  the  equilibrium  constant.  These  data  are  summarized  for  the  two 
isotopes  in  Fig.  5  for  both  annealed  and  cold  rolled  single 
crystals.  Trapping  in  contrast  to  rapid  diffusion  along 
dislocation  cores  is  consistent  with  these  results. 

IV.  THEORETICAL  STUDY  OF  THE  TRAPPING  OF  HYDROGEN  BY 

DISLOCATIONS 

Min  Yan,  a  graduate  student  at  Michigan  Tech,  applied  a 

28 

modification  of  the  Embedded  Atom  model  (EAM  )to  calculate  the 
binding  of  a  hydrogen  atom  in  an  interstitial  site  of  a-iron.  In 
general,  these  workers  at  Sandia  National  Lab  had  not  carried  out 
any  calculation  on  BCC  metals.  Although  the  energy  value  for 

Xiaolin  Shi,  Thesis  at  Michigan  Tech,  1989 

28 

M.  Baskes,  C.  F.  Melius  and  W.  D.  Wilson,  Zeit.  f\ir  Phys. 

Chem. ,  Neue  Folge,  116  (1979)  289.  See  also  Proc.  Effect  of 

Hydrogen  on  the  Behavior  of  Materials,  Jackson  Lake  Lodge, 

Wy. 1980 


14 


THE  DISLOCATION  CONFIGURATION 


Fig.  6 


Placed  perpendicular  to  the  Dislocation  line.  Cross  marks  the  approximate 
location  of  the  mimimum  energy  of  a  Hydrogen  Atom. 


•  NH* 


Perspective  View 


rt^AX  •  3.f 


P(o».  CeMeuf  plot  •t  Ky^o9#n  \f>  rp^\or% 

•  hOvfl  %f%  r»f.  4<b». 


miN  •  -4.T 


Fig.  7  Two  Gsmplementary  Views  of  the  Variation  of  the  Potential  of  A 

Hydrogen  Near  the  Core  of  A  Dislocation  in  Iron  Over  the  XY  Plane 
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binding  at  an  interstitial  site  was  quite  high  it  agreed  well  with 
a  calculation  made  by  another  method. 

The  next  step  she  took  was  to  introduce  a  (100)  edge 
dislocation  into  the  lattice  by  movement  of  half  planes  with  a 
common  slip  plane.  The  relaxation  of  the  atoms  in  the  vicinity  of 
the  core  was  done  initially  using  the  elastic  strains  of  a 
continuum  model  of  a  dislocation.  Then  displacements  due  to  the  x 
and  y  contributions  to  the  interatomic  forces  were  calculated  from 
the  modified  EAM.  All  the  atoms  were  treated  pairwise  through 
several  iterations.  The  region  which  was  considered,  extended  to 
five  lattice  parameters  from  the  core. 

The  displacements  are  Indicated  in  Fig.  6  and  the 

bending  of  the  slip  plane  Is  evident.  Isopotential  energy  lines 

for  a  hydrogen  atom  in  the  vicinity  of  the  core  are  shown  in  two 

different  versions  in  Fig.  7.  The  minimum  for  hydrogen  is  located 

Just  below  the  slip  plane.  The  influence  of  the  iron  atom  just 

behind  the  XY  plane  of  the  contour  plot  is  not  clearly  revealed. 

The  variation  of  the  potential  along  the  dislocation  core  within 

the  additional  half  plane  of  the  dislocation  is  shown  in  the 

Fig.  8.  The  indicated  binding  energy  is  0.35  eV  which  is  larger 

29 

than  the  value  Park  deduced  in  his  thesis  ,  namely  37  kJoules 

for  edge  dislocations  and  26  for  screw  dislocations.  The  latter 

30 

value  is  in  good  agreement  with  the  value  Hirth  reported,  namely 
20-30  kJoules  for  screw  dislocations.  In  the  right  hand  panel,  the 
potential  is  plotted  for  jog  created  by  removing  two  atoms  from 
the  half  plane.  As  noted  in  the  iegend,  the  calculated  potential 
well  is  about  0.05  eV  deeper  when  the  jog  is  included. 

Incomplete  Study  of  Dislocations  in  FCC  Metals 

Eric  Bylaska,  an  undergraduate,  undertook  a  study  of 
equivalent  calculation  for  hydrogen  in  FCC  palladium  after 
completing  an  undergraduate  thesis  in  Physics.  However,  it  took 
much  more  coding  to  obtain  a  block  containing  a  (111)  [110] 


29 

30 


See  Ref.  3. 
See  Ref.  8. 
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dislocation  and  with  a  sufficient  number  of  atoms  than  with  the 
(100)  [010]  dislocation  of  iron.  The  direction  of  viewing  the  FCC 
palladium  is  along  a  <112>  direction.  The  movement  of  the  half 
planes  to  create  an  edge  dislocation  is  indicated  in  the  three 
portions  of  Fig.  10.  Unfortunately  the  elastic  displacements  have 
not  yet  been  computed. 

POSITRON  TRAPPING  IN  POLYMERS 
31 

Ronghal  Wang  studied  the  trapping  of  positrons  in  a  series 
of  epoxy  resins.  The  Initial  purpose  was  to  assess  the  thermal 
variation  of  the  free  volume  of  the  polymers  and  the  glass 
transition  temperature.  In  addition,  the  effect  of  changes  in  the 
structure  of  the  polymers  on  the  free  volume  was  also  studied; 
specifically,  changes  in  the  fvmctionality  of  the  amino  group  and 
changes  in  the  number  of  CH^  in  the  side  chains  of  the  curing 
agents  were  made.  N-methylation  of  the  amino  group  is  designated 
C2M0,  C2M1,  and  C2M2;  these  symbols  mean  that  0,  1  and  2  methyls 
have  been  added  to  diamino-ethane.  Thus  the  number  of  functional 
N-H  bonds  which  can  react  with  the  epoxy  groups  was  reduced.  The 
symbols  C3M0  and  C4M0  signify  that  n-propyl  diamine  and  n-butyl 
diamine  were  used  as  curing  agents. 

It  was  possible  to  measure  the  coefficient  of  expansion  of 
the  free  volume  over  a  range  of  temperatures.  In  the  literature, 
it  is  postulated  that  the  free  volume  is  constant  and  possibly  is 
even  negligibly  small,  just  below  the  glass  transition  temperature 
T  .  Using  the  lifetime  of  the  so-called  third  or  longest 
component  of  the  lifetime  spectrum  t^,  it  was  found  that  the 
fractional  free  volume  varied  from  3.7  to  5.4  per  cent  at  T  for 

9 

the  five  polymers  and  decreased  to  2  to  3  per  cent  at  -50°C. 
Contrary  to  the  published  assumptions,  the  free  volume  was  not 
constant.  When  the  data  were  extrapolated  to  absolute  zero,  the 
residual  volume  was  -0.34  ±  0.58  percent.  This  indicates  the 
inherent  problem  of  extrapolating  the  data  over  a  long  range, 
i.e.,  more  than  200°. 


Rong-Hai  Wang,  Thesis  Michigan  Tech  1991 
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Fig.  lOa  Fcx:  Lattice  Viewed  Along  [lT2)  Direction 


Fig.  10b  Same,  with  two  adjacent  <110>  planes  removed 


The  free  volume  did  increase  significantly  above  T  .  While 
other  researchers  had  used  such  a  change  in  the  slope  of  the 
positron  lifetime  curve  to  find  T  ,  no  comparison  with  the  thermal 

9 

expansion  of  the  solid  as  a  whole  was  made  nor  was  the  free  volume 
measured  at  such  low  temperatures.  In  addition,  the  thermal 
dilation  of  the  occupied  volume,  i.e.,  the  C-H  chains  themselves, 
was  estimated  by  the  increase  in  the  second  component  Both  of 
these  lifetimes  were  converted  into  fractional  free  volumes  from 
the  products  t  I  and  t  I  (where  I  is  the  intensity  or  better, 

2  2  3  3 

the  fraction  of  the  positrons  annihilating  in  a  specific  manner) 

32 

using  an  expression  Nakanishi  and  Jeans  had  developed.  As  shown 
in  the  table  below,  the  sum  of  the  two  volume  increases  was  in 
good  agreement  with  the  thermal  expansion  coefficients  measured  in 
a  torsional  apparatus.  Because  of  the  possibility  of  the  stylus 
distorting  the  heated  sample,  the  value  in  the  last  column  tends 
to  be  high,  except  for  C2M2  which  is  only  one  of  the  five  polymers 
to  melt  well  above  T  .  It  is  interesting  that  most  of  the  dilation 

9 

below  T  ,  in  our  experiments,  comes  from  the  free  volume  and  not 

9 

from  the  occupied  volume  as  previously  assumed. 

In  view  of  this  inaccuracy  in  measuring  the  dilation 
mechanically,  there  is  good  agreement  between  the  measured  values 
and  those  deduced  from  the  positron  lifetimes.  One  conclusion 
drawn  is  that  the  nature  of  the  "joint"  between  the  curing  agent 
segments  and  the  epoxy-molecule  plays  a  minor  role  compared  with 
the  composition  of  the  resin,  i.e.,  the  mixing  ratio  needed  for 
complete  reaction. 


H.  Nakanishi  and  Y.  C.  Jeans,  in  Positrons 
and  Positronium  Chemistry,  (1988)  pp.  159-162. 
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TABLE  1.  COMPARISON  OF  THE  THERMAL  EXPANSION  BEHAVIOR 


OF  THE  FIVE  EPOXY-RESINS 


Below  the  Transition  Temperature 


Symbol 

Sum  of  2  Alphas 

Exper.  Dilation 

C2M0 

16.35  X  10'® 

17.5  X  10'® 

C2M1 

16.51 

21.5 

C2M2 

17.53 

20.1 

C3M0 

18.19 

18.3 

C4M0 

14.01 

20. 1 

Above  Transition 

Temperature 

C2M0 

50.45 

61.9 

C2M1 

59.0 

55.3 

C2M2 

66.5 

40.8 

C3M0 

59.5 

40.1 

C4M0 

42.1 

40.8 
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GOAL  AT  THE  TIME  OF  FUNDING: 

Use  positrons  to  sample  whether  Hydrogen  traps  were 

OCCUPIED  AND  HOPEFULLY  TO  DETERMINE  THEIR  FRACTIONAL 
OCCUPATION. 


In  addition  to  this  being  successfully  ACHIEVED, 
THE  FOLLOWING  RESULTS  WERE  OBTAINED: 

1.  A  NON-DESTRUCTIVE  METHOD  OF  DETERMINING  THE 
CONCENTRATION  OF  EDGE  AND  SCREW  DISLOCATIONS  WAS 
DEVELOPED,  AS  DESCRIBED  ABOVE, 

2.  A  THEORY  OF  THE  THERMAL  REMOVAL  OF  DISLOCATIONS 

WAS  VERIFIED.  J.  C.  M.  Li’S  MECHANISM  INVOLVES  THE 
ANNIHILATION  OF  DISLOCATION-DIPOLES  BY  GLIDE.  THIS  IS 
IN  CONTRAST  TO  THE  PREVIOUSLY  ACCEPTED 

VACANCY-DRIVEN  DIFFUSION  MECHANISM,  HOWEVER,  SINCE  THE 
VACANCY  CONCENTRATION  COULD  BE  INDEPENDENTLY  DETERMINED 
BY  THE  POSITRONS,  IT  WAS  CLEAR  THAT  THERE  WERE 
INSUFFICIENT  NUMBERS  OF  VACANCIES  TO  DIFFUSE  TO  THE  CORE 
OF  DISLOCATIONS  AND  REMOVE  THEM  BY  CLIMB. 

3.  Mechanism  of  dislocation  generation  yields  a 

BALANCED  NUMBER  OF  POSITIVE  AND  NEGATIVE  DISLOCATIONS 
TO  FORM  DIPOLES  ON  THE  SAME  PLANE  '  THAT  IS,  THE  EXCESS 
OF  ONE  SIGN  OVER  THE  OTHER  IS  ONLY  A  FEW  PERCENT  OF  THE 

TOTAL.  This  was  also  deduced  from  the  annealing 
STUDIES  using  Ll'S  THEORY. 

4.  A  method  was  developed  for  extracting  isothermal 
KINETIC  analysis  INFORMATION  FROM  ISOCHRONAL  ANNEALING 
RUNS.  A  COMPARISON  SHOWED  GOOD  AGREEMENT  BETWEEN  THE 
TWO  METHODS.  MaINLY  USING  ISOTHERMAL  RUNS,  IT  WAS 
ESTABLISHED  THAT  THE  RATE  OF  DISLOCATION  ANNIHILATION 
WAS  A  SECOND  ORDER  REACTION  AND  DEPENDED  ON  THE 
CONCENTRATION  OF  BOTH  OF  THE  DIPOLE  COMPONENTS. 
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5.  The  activation  energy  for  screw  dipole  migration 

WAS  SMALLER  THAN  THAT  FOR  EDGE  DIPOLE  MIGRATION  -  AND 
BOTH  WERE  SMALLER  THAN  178  K JOULES  FOR  VACANCY  MIGRATION 
-SELF-DIFFUSION  -  IN  ALPHA-IRON, 

5.  Dislocation  profile  formed  by  mechanical  polishing 
WAS  established  -  Concentration  is  a  decreasing 

LOGARITHMIC  FUNCTION  OF  DEPTH.  RESULTS  ARE  CONSISTENT 
WITH  THE  STUDY  OF  DISLOCATION  REMOVAL  DURING  CREEP 
IN  BOTH  BCC  AND  FCC  METALS  BY  ArSENAULT  AND  COWORKERS. 

[  TWO  ACTIVATION  ENERGIES  AND  A  DECREASING  PROFILE  1 

6.  Direct  influence  of  dislocations  in  the  surface 

LAYERS  ON  BULK  PERMEATION  OF  HYDROGEN  IN  IRON  WAS 
MEASURED.  WHEN  COUPLED  WITH  EFFECTS  OF  MECHANICAL 
POLISHING  OF  DIFFUSION  SAMPLES,  THIS  OBSERVATION  OFFERS 
A  PARTIAL  EXPLANATION  OF  THE  LARGE  VARIABILITY  OF 
DIFFUSIVITY  MEASUREMENTS  BELOW  300°  C.  ABOVE  300  ALMOST 
OF  THE  EDGE  AND  SCREW  DISLOCATION  ANNEAL  OUT  DURING  THE 
EXPERIMENT  AND  HENCE  HAVE  LITTLE  EFFECT  DN  DIFFUSION.. 

7.  Measurement  of  the  diffusivity  of  deuterium  and 

HYDROGEN  INDICATE  A  NON-CLASSICAL  VALUE  OF  THE  EFFECT  OF 
ISOTOPIC  SUBSTITUTION, 

8.  Theoretical  computation  of  the  diffusion  activation 

ENERGY  AND  ENERGY  CONTOURS  FOR  HYDROGEN  IN  A 
DISLOCATION  AND  A  JOG  OF  BCC  IRON 

9.  Positron  lifetime  measurements  on  a  series  of  epoxy 

RESINS  CAST  DOUBTS  ABOUT  THE  PUBLISHED  ASSUMPTIONS  ABOUT 
THE  BEHAVIOR  OF  THE  FREE  VOLUME  BELOW  THE  GLASS  TRANSITION 
TEMPERATURE.  IT  IS  NEITHER  SMALL  NOR  CONSTANT. 
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HYDROGEN  PERMEATION  and  POSITRON  ANNIHILATION 


STUDY  of  ALPHA  IRON  CRYSTALS* 
by 

James  T.  Waber 

Physics  Department 
Michigan  Technological  University 
Houghton,  MI  49931-1295 

Dislocation  densities  in  deformed  iron  single  crystals  were  determined  by  positron 
annihilation  techniques.  Electrochemical  measurements  of  hydrogen  diffusion  were  made 
on  the  same  crystals.  It  was  shown  that  apparent  diffusivity  D^  was  inversely  proportional 
to  the  dislocation  density  -  results  cover  three  different  modes  of  deformation.  Mechanical 
polishing  of  annealed  specimens  with  even  fine  abrasives  introduced  sufficient  dislocations 
to  reduce  the  diffusivity  tenfold. 

Cathodic  deposition  of  hydrogen  and  of  deuterium  atoms  followed  by  rapid  removal 
by  potential  reversal  was  studied  on  similar  crystals.  Holdup  of  the  hydrogen  atoms  in 
surface  layers  was  increased  by  the  dislocations.  The  effect  was  stronger  with  deuterium 
charging.  Surface  coverage  was  approximately  5  percent  of  the  available  sites. 
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To  be  presented  at  the  UM-ARO  Workshop  on  Surface  Science 
and  Technology,  Ann  Arbor,  MI  on  November  7-9,  1990. 


Hartree-Fock  Calculations 


on  N-Alkanes  Adsorbed  on  Graphite* 
by 

James  T.  Waber 
and 

Max  Seel 

Department  of  Physics 
Michigan  Technological  University 
Houghton,  MI  49931-1295 

Scanning  tunneling  microscope  study  by  McGonigal,  jl.  of  simple  aliphatic 
hydrocarbons  on  pyrolytic  graphite  shows  highly  ordered  2-D  pattern  of  the  adsorbed 
molecules.  The  zig-zag  C-C  chain  is  apparently  slightly  elongated  to  be  commensurate  with 
graphite  cell  spacing. 

Unrestricted  Hartree-Fock  calculations  were  carried  out  with  an  optimized  basis  set 
of  Huzinaga^  for  a  graphite  cluster  as  well  as  for  several  small  alkanes.  This  set  had  been 
used  to  determine  correctly  the  crystal  structure  of  several  normal  alkanes  and 

change  in  crystal  symmetry  which  occurs  when  n  is  either  even  or  odd.^ 

The  effect  of  superposing  the  alkane  on  the  graphite  surface  was  studied  and  the 
location  varied  to  find  a  local  energy  minimum. 

Reasons  how  the  adsorbate  enhances  the  tunneling  current  from  specific  graphite 
sites  will  be  discussed. 

1.  G.C.  McGonigal,  R.H.  Bernhardt,  R.H.  and  DJ.  Thompson,  App.  Phys.  Lett,  (1) 
1990,  28-30. 

2.  S.  Huzinaga,  "Gaussian  Basis  Sites  for  Molecular  Calculation,"  Elsevier,  Amsterdam- 
New  York  (1984). 

3.  D.J.  Wyble,  M.  Seel  and  J.T.  Waber,  Solid  State  Comm.,  M  (5)  (1987)  827-829. 
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To  be  presented  at  loth  International  Conference  on 
Scanning  Tunneling  Microscopy,  Interlaken  Switzerland, 
August,  1991. 


FREE  VOLUME  CHANGES  AND  POSITRON 


LIFETIMES  IN  EPOXY-DIAMINE  POLYMERS 
by 

Ronghai  Wang  and  James  T.  Waber 
Physics  Department 
and 

John  Williams 

Chemical  Engineering  Department 


MICHIGAN  TECHNOLOGICAL  UNIVERSITY 
Houghton,  Michigan 


Abstract 


Five  epoxy-cured  diamine  resins  have  been  studied  with  positron 
annihilation  techniques.  The  data  were  fitted  by  three  lifetimes. 

Principal  attention  was  paid  to  the  third  and  longer  lifetime  which  is 
associated  with  ortho-positron ium  (0-  Ps) .  Temperatures  ranged  from  -50 
to  150*  to  C.  Changes  in  the  volume  -  coefficient  of  expansion  and  the 
glass-transition  temperature  were  determined.  It  was  established  that 
aging  at  temperatures  for  periods  several  times  longer  than  a  measurement 
required  had  no  discernible  effect  on  these  parameters. 

The  correlation  between  the  temperature  variation  of  the  free  volume  and 
that  derived  from  WLF  theory  of  viscosity  was  investigated. 


To  be  presented  at 
First  Pacific  Rim  Conference  on 
Advanced  Materials  and  Processing 


Hangzhou  China 
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